A black hole solution to the cosmological monopole problem 
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We propose a solution to the cosmological monopole problem: Primordial black holes, produced in 
the early universe, can accrete magnetic monopoles before the relics dominate the energy density of 
the universe. These small black holes quickly evaporate and thereby convert most of the monopole 
energy density into radiation. We estimate the range of parameters for which this solution is possible: 
under very conservative assumptions we find that the black hole mass must be less than lO^gm. 
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Magnetic monopoles are a generic prediction of Grand 
Unified Theories (GUTs). Monopoles always appear as 
stable topological entities in any GUT that breaks down 
to electromagnetism [1,2]. The monopole mass is M^, ~ 
ry/e, where rj is the energy scale of the phase transition 
and l/e — y/137. Hence, for r] ~ 10^^ GeV, monopoles 
of mass Aim — 10^^ GeV should exist. We use this mass 
scale as our fiducial value throughout, for the sake of def- 
initeness, but note that an analogous calculation can be 
done for any other scale*. The magnetic charge of the 
monopoles, g, is determined by the quantization condi- 
tion eg = n/2 {n is an integer). The abundance of these 
monopoles is an open question. The Kibble mechanism 
[3] predicts roughly one monopole per horizon volume at 
the time of the GUT phase transition (due to causality 
constraints). However, this estimate provides a severe 
overabundance of the number of monopoles: they over- 
close the Universe by many orders of magnitude. Since 
they are stable for topological reasons, they do not de- 
cay; it has also been shown that the rate of pair annihila- 
tion cannot reduce the monopole number sufficiently [4] . 
Their abundance is expected to be too high to be consis- 
tent with astrophysical observations of magnetic fields in 
our Galaxy [5] and consequences of catalyzed nucleon de- 
cay in stars [6,7]. This overabundance of GUT monopoles 
is known as the monopole problem. 

An inflationary epoch [8] may reduce the density of 
monopoles in the Universe (for some alternative ap- 
proaches see e.g. [9]). The rapid expansion phase di- 
lutes the monopole abundance. The re-heating tempera- 
ture must be sufficiently low in order not to produce the 
monopoles again. On the other hand, since any initial 
baryon-antibaryon asymmetry will be also washed out, 
the re-heating temperature must be sufficiently high in 
order to provide a mechanism for the observed baryon- 
antibaryon asymmetry. Here we explore new solutions to 
the monopole problem which, for example, would obviate 



*For example, for rj ~ lO^^GeV (which is favored by su- 
persymmetric extensions of the standard model and by ex- 
perimental constraints on the proton lifetime), the monopole 
mass would be 10^* GeV. 



the need for a low re-heating temperature. 

As the universe expands, the cosmological horizon 
grows to encompass many monopoles. The predicted ini- 
tial ratio [10] between the monopole and entropy density 
is given by 
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where p ~ 0.1 is a geometrical factor and ^ 100 is 
the number of effectively massless degrees of freedom 
at a given temperature. This ratio would remain con- 
stant as the universe adiabatically expands. However 
non-adiabatic processes, such as monopole/antimonopole 
annihilation or accretion by black holes, reduce this ratio. 

Here we first review the relic monopole abundance sub- 
sequent to annihilation effects only, and later include a 
new effect: the fact that primordial black holes can eat 
monopoles. The relic abundance is far too large with 
annihilation only, but this problem can be resolved with 
black holes. 

The motion of a monopole in a plasma is basically a 
random walk where the typical step is the mean free path. 
During each step a monopole moves with thermal velocity 
{2T/Mmf''^. The mean free path / depends on the 
scattering of a monopole on charged particles in plasma 
and is given by [10]: 
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The factor h depends on the number of degrees of free- 
dom of light charged particles and its value is 6 ~ 100. 
Monopole-antimonopole capture occurs when g^/r ~ T, 
which gives the capture radius Tc ~ 9^ /T. Annihilation 
effectively takes place as long as the mean free path of a 
monopole is smaller than the capture radius for annihi- 
lation [11]. This implies that the annihilation practically 
stops at the final temperature 
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10"GeV, 
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and the monopole to entropy density ratio freezes at the 
final value [4,10] 
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The annihilation mechanism is effective only if the ini- 
tial monopole to entropy density ratio is greater than 
the one given by (4). Otherwise, the initial ratio remains 
constant as the universe expands. For our fiducial choice 
of parameters with ~ 10^® GeV, annihilation is ir- 
relevant and the initial monopole to entropy density ra- 
tio of Eq.(l) remains constant; for other parameters and 
monopole masses, however, Eq.(4) should be used. The 
entropy density at a given temperature T is s ^ g^,T^. 
The total mass of monopoles within the Hubble volume 
is thus 

(5) 

where Vh ~ is a Hubble volume at a given 

Hubble time th = Q.2,g7^''^ Mpi/T"^ . This mass will 
equal the total mass within the Hubble volume [Mh ~ 
0.357^^^M|,,/T2) for T ~ lO^GeV. Thus, monopoles 
dominate the energy density of the Universe at T < 
lO^GeV. Thus the predicted present day monopole flux 
grossly violate many late-time cosmological constraints, 
if one considers annihilation as the only mechanism for 
reducing the monopole number. 

Now we will consider the effects of black holes in re- 
ducing the monopole number. The early universe can 
produce large numbers of primordial black holes, via a 
number of processes [12-15]. The earliest mechanism for 
black hole production can be fluctuations in the space- 
time metric at the Planck epoch. Large number of pri- 
mordial black holes can also be produced by nonlinear 
density fluctuations due to oscillations of some (scalar) 
field. If within some region of space density fluctuations 
are large, so that the gravitational force overcomes the 
pressure, we can expect the whole region to collapse and 
form a black hole. Black holes can also be produced at 
first and second order phase transitions in the early uni- 
verse [14]. Collapse of cosmic string loops and closed 
domain walls also yields black holes. 

We will show that black holes could solve the monopole 
problem by capturing them before they become danger- 
ous for the standard cosmology. Since the black holes 
evaporate, much of the dangerous excess monopole en- 
ergy density would then be converted into harmless ra- 
diation. If the universe contains enough black holes to 
remove monopoles quickly enough, but not so many of 
them to caiise deviations from standard cosmology, then 
the monopole problem would be solved. 

Primordial black holes can have a variety of masses. A 
typical mass would be that of the mass inside the horizon 
at the time of formation, so that the mass would range 
roughly from Mpi (black holes formed at the Planck 
epoch) to Msun (black holes formed at the QCD phase 
transition). However, the mass could be much smaller 
[13] . For example, black holes formed by collapsing closed 



domain walls and string loops are typically much lighter 
than a horizon mass. After formation black holes can 
evolve further by evaporating, merging with each other 
and accreting surrounding matter. Since the number of 
black holes of a given mass within a horizon volume as 
a function of time is strongly model dependent, we do 
not limit ourselves to any particular model. We instead 
solve the dynamical equation for the time evolution of 
the number density of monopoles with free parameters 
that determine the abundance of black holes, and then 
estimate the range of these parameters that is required 
for the solution of the monopole problem. We will show 
that black holes of mass M^h < lO^'gm (a mass range 
that does not usually put any serious constraints on stan- 
dard cosmology) are dynamically capable of removing the 
monopoles. 

Monopoles move with non-relativistic velocities 
vm << c in the early universe, since they are magneti- 
cally charged and their velocity is damped due to inter- 
action with surrounding plasma. The cross section for 
gravitational capture of a non-relativistic massive par- 
ticle (monopole) by a black hole of gravitational radius 
Rbh is [16]: 

C7g = inic/vM^Rlh- (6) 

The monopole flux is given by rimVm^ where iim is the 
number density. If there is only one black hole in the 
Hubble volume, the number of monopoles that get cap- 
tured by the black hole, during the Hubble time r/j, is 
nmVm<^gTh- The total number of captured monopoles is 
obtained by multiplying this quantity by the number of 
black holes within the Hubble volume, ribhVh, where Vh 
is the Hubble volume. This number should be compared 
with the total number of monopoles contained within the 
Hubble volume, nmVh, thus giving the condition 

nbh'JgVmTh > 1 , (7) 

This inequality represents the condition under which the 
characteristic time for monopole capture is shorter than 
the Hubble time. We will now ascertain the resulting 
number density of monopoles by examining the dynamics 
in more detail. As we will see, the term on the left hand 
side of (7) will play the role of the suppression factor 
for the evolution of the monopole number density. How 
big this factor should be in order to solve the monopole 
problem, we will see from the equation that describes the 
dynamics of the process. 

The dynamical equation (i.e. the detailed balance 
equation) for the time evolution of number density of 
monopoles in comoving volume is 

= -^"m - nbhf^gVmrim - 3-n„ . (8) 

The first term on the right hand side is present as 
long as the monopole-antimonopole annihilation is tak- 
ing place. The parameter A was estimated in [11] as 
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A = g^/ibT^) = g'^gV'^t/{QMMpi) w m/Mpi. The sec- 
ond term takes into account the monopole capture by the 
black holes while the third term is due to expansion of 
the universe (a is the scale factor of the universe) . 

We now examine the term richer gVm- The energy den- 
sity of the black holes at any given time must be less 
than the energy density of the universe, so that pf,h = 
MbhUbh ~ f9*T'^, where / < 1, while M^h is the mass 
of a typical black hole at that scale. The number den- 
sity of black holes is thus ribh ^ fg^T'^/Mbh- For a given 
black hole mass M^h, the radius Rbh = 2Mi,?i/M|>(. The 
relative black hole-monopole velocity is just a monopole 
random walk velocity = vt/VN, where vt is a ther- 
mal velocity while N is the number of random walk steps. 
At some distance r, we have \/]V = r/l, where I is the 
length of the step (a mean free path). The characteris- 
tic distance is the gravitational capture radius calculated 
from the condition MbhMm/{Mpirc) ~ T. This gives 
VM = V2MpJ{bMbhMm)- The mass of a typical black 
hole is just a fraction of the mass within the horizon at 
any given time. In other words, 
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where /3 < 1. In the absence of complete knowledge of 
the distribution of black holes as a function of time, an 
assiimption must be made about the time dependence of 
/ and p. In our analysis of the dynamics, we will assume 
that the parameters / and (3 do not change significantly 
with time as they describe some average features of the 
system^. Substituting these relations into UbhO'gVm we 
obtain 

nbhTgVm = fP'^bMm « O.lfp^Mpi . (10) 
The dynamical equation (8) now reads 

drim lot 2 nwfl2,^ 3 

-dT = " M^"™ ~ - 2^"" • ^^^^ 
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The solution to this equation is 

g-0.1//3=t/tpi 



Cii3/2 + (0.l7r/M3//32)i/2t3/2erf y^0.1f(3H 
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where erf(x) = 2/y/Tr e~*^ dx and Ci is a constant of 
integration. 

The gravitational capture of monopoles by the black 
holes effectively takes place as long as the capture radius 



^We note that this assumption requires black holes to form 
at a range of temperatures or to grow due to accretion. 



is larger than the mean free path of monopoles. This 
condition is satisfied for T < « 
{Q.lpi^Y/^Mpi. While the electromagnetic monopole- 
antimonopolc capture stops at some finite temperature, 
the gravitational capture starts at some high tempera- 
ture and practically never ends. This is the consequence 
of the fact that the universe can accommodate larger and 
larger black holes as time passes, so that the gravitational 
capture radius grows faster than the mean free path of 
monopoles. Hence wc can neglect the monopole annihi- 
lation contribution in eq. (12) (the second term in de- 
nominator). We note that in general dropping the second 
term is a conservative assumption in the sense that we 
are underestimating the destruction of monopoles. Thus 
we have 
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where the initial ratio between the monopole and entropy 
density is given by cq. (1). We sec that the gravitational 
capture of the monopoles by the black holes exponentially 
reduces this ratio. We note that the most massive black 
holes contribute the most, while the contribution of much 
lighter ones can be neglected. Also, if ti <^ t we can 
neglect the term ti in the exponential. This will be true 
in practically all the cases of interest here. 

Now, we can derive the values of / and /3 needed for a 
solution of the monopole problem. The first bound on the 
monopole abundance comes from the requirement that 
the universe must be radiation dominated at nucleosyn- 
thesis, i.e. T ~MeV. This implies that the ratio between 
the monopole and entropy density can be at most 
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We can derive a more stringent bound by requiring 
that the mass density of the monopoles today should not 
exceed the critical density: Mmfim < Pc- This yields 
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where we used the current value of the entropy density 
So = lO'^cm"'^. There are a number of other late time 
constraints [5]. For example, since monopoles are mag- 
netically charged, they can be accelerated by galactic 
magnetic fields. In this process the galactic magnetic 
field gets dissipated. In order to avoid a conflict with ob- 
servations, the upper limit on monopole number density 
is 



rim < 10" 
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Certainly, the strongest constraint comes from the con- 
sequences of monopole induced baryon catalysis in neu- 
tron stars [7]: 
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This implies the bound on a current ratio between the 
monopole and entropy density: 
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Substituting eq. (1) and eq. (18) into eq. (13) we can 
derive 



fpHf/tpi > 100, 



or, equivalently. 



f(3^ {m,pi/Tff > 100 



(19) 



(20) 



for our fiducial value of Mm — 10^^ GeV. This is the 
necessary condition for a solution of the monopole prob- 
lem. Here, tf and T/ denote time and temperature re- 
spectively when the process of gravitational capture of 
monopoles by black holes stops. We now discuss whether 
or not the universe can reasonably produce black holes 
that satisfy this condition. 

Black holes more massive than 10^^ gm have lifetimes 
longer than the present age of the universe. While these 
black holes are dynamically capable of swallowing the 
monopoles, the energy density that was in the monopoles 
remains in the form of mass rather than radiation and 
hence overdoses the universe. Hence we must consider 
lighter black holes. 

Black holes of mass Mbh < lO^^gm evaporate before 
the present epoch, thereby converting the monopole en- 
ergy density to radiation which redshifts away. However, 
observational constraints from earlier epochs have been 
previously examined by other authors [15] and found to 
be more severe. These constraints are sensitive to the 
assumptions in the model, e.g., model of formation of 
the black holes, equations of state in different epochs 
of the early universe, whether or not the universe expe- 
riences a period (or multiple periods) of inflation, etc. 
Ref. [15] shows that black holes with mass over lO^g 
could not have been produced with appreciable abun- 
dance relative to the total energy density without vio- 
lating observational constraints due to the cosmic mi- 
crowave background (CMB), nucleosynthesis, geometry 
of the universe, etc. On the other hand, black holes less 
massive than lO^gm are not seriously constrained by ob- 
servations: they evaporate completely within G.lsec, i.e. 
before nucleosynthesis takes place, and are compatible 
with results of CMB experiments. 

Thus, we work with the constraint that only black 
holes of Mhh < lO^gni can be responsible for removing 
monopoles from the universe. Using Eqs.(20) and (9), 
together with this lower bound on the black hole mass, 
we find 
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As an example, we will consider Tf — lO^GeV, or, 
equivalently tf ^ lO^^ip;, for the final temperature. 
Then the condition (19) becomes 



We plot the allowed range in Fig. 1. 
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FIG. 1. The allowed parameter range for the solution of 
the monopole problem (for Mm = lO^^GeV) lies above the 
curve 

/?= (10-16//)'''^ (shaded region on the Log- Log plot). 
Parameter / is the fraction of the total energy density of the 
universe that must be in black holes, while /3 is the fraction 
of the horizon mass that represents a typical (average) black 
hole at that time. Even for values of / <C 1 we can have 
reasonable values for (3. 

For example, if / ~ 10^^ (1% of the energy density of 
the universe is in black holes), the allowed range for f3 is 
(3 > 10~^, which means that the mass of a typical black 
hole at that time should be at least 10"^ of the horizon 
mass. This is not difficult to envision. The horizon mass 
at T ~ lO^GeV is roughly Mh ^ lO^^gm, so the average 
black holes mass should be M^h > lO^gm. Note that 
the number of black holes within a horizon is roughly 
given by Ni,fi = ///3, but since /3 gives only an average 
description of the system, the real number of black holes 
could be quite different from Nbh defined in this way. We 
therefore do not use N^h as a reliable quantifier. 

We have shown that, as far as the dynamics of the 
process is concerned, a reasonable number of primordial 
black holes is capable of capturing monopoles from the 
early universe. We now discuss problems that might po- 
tentially arise. 

After absorbing monopoles, a black hole will be- 
come charged. In our calculation we assumed that the 
monopole capture by the black hole is charge blind, i.e. 
we neglected the electromagnetic interaction between the 
monopole and black hole. Now, we justify this assump- 
tion. Let's consider a black hole which has already cap- 
tured Nm monopoles and antimonopoles. Most of the 
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monopolcs and antimonopolcs inside the black hole will 
annihilate, leaving monopoles inside (due to mag- 

netic charge fluctuations). The gravitational mass of the 
black hole is enhanced to M^h + NmMm] the magnetic 
charge of the black hole becomes \/Nmg, where 5 is a 
magnetic charge of the monopole. The gravitational force 
between this black hole and a monopole outside it is then 

Fg = G{Mbh + NmMm)Mm/r^ . (23) 

The corresponding magnetic force is 

Fm = ^N^g^/r^ ■ (24) 

If Nm is large {Nm > 10^^), gravitational force is always 
greater than the magnetic. Even if Nm is small, one can 
easily check that the gravitational attraction between a 
monopole of mass lO^^GeV (with a unit charge g) and 
a black hole that captured Nm monopoles is larger than 
the electromagnetic repulsion between them for any black 
hole with 

Mbh > 10^^/N^Mpl, (25) 

a condition which is easily satisfied for the black holes of 
interest. 

Another problem that might arise in this scenario 
can come from magnetic charge fluctuations: a black 
hole that captures Nm monopoles and anti-monopoles 
will typically have accumulated a residual net magnetic 
charge of s/N^. The question is, what happens to 
this magnetic charge? One possibility would be that 
black holes Hawking-radiatc monopoles. A monopole is a 
highly coherent state of many gauge quanta and emission 
of a monopole by a black hole will be highly suppressed. 
Even if this process is somehow allowed [17], the Hawk- 
ing temperature of a black hole becomes of order of a 
monopole mass only at tlic^ end of evaporation, and a 
black hole could radiate much less monopolcs than orig- 
inally eaten. Gcnerically, since the Hawking radiation 
can not violate the gauge symmetry, a black hole can 
not evaporate completely. Instead it leaves a remnant 
- an extreme magnetically charged Reissner-Nordstrom 
black hole. The mass of the (non-rotating) remnant Mr 
must be greater than the magnetic charge Qm of the 
black holes, Mj. > Qm (in appropriate units), or oth- 
erwise the remnant would be a naked singularity. One 
has to check whether these remnants violate some of the 
observational constraints. Due to the attractive force be- 
tween the remnants with opposite charges, their number 
density will be reduced (charges annihilate and the rest of 
the mass evaporates freely). However, here we adopt the 
most conservative scenario in which all the black holes 
which participate in accreting monopoles leave massive 
remnants. The mass density \/Nm that remains trapped 
in the extremal Reissner-Nordstrom black holes could po- 
tentially be problematic. The energy density in these 



remnants must be less than that of ordinary radiation at 
the time of primordial nucleosynthesis, in order for pre- 
dictions of element abundances to be unaff'ected. Here 
we will compare the residual number of monopoles inside 
black holes and ensure that it is smaller than the num- 
ber of monopoles allowed by the nucleosynthesis bound 
in Eq.(14). The total number of monopolcs within the 
horizon at temperature T is Ntot = Mtot/Mm, or, using 
Eq.(5), 

A^.ot~10-^(^)'. (26) 

We note that this number is proportional to rim/s. If 
all of the monopoles within horizon are eaten by black 
holes at some temperature (e.q. T ~ lO^GeV), then 
the remnant monopole mass inside the black holes is 

\/Nm oc ^ (^)j given in Eq.(l). On the other hand, 
the allowed number of monopoles is proportional to 
A^aiiowed OC (nm/s)MeV given in Eq.(14). Thus, we can 
calculate that at T ~ lO^GeV, there are ~ 10^^ 
monopoles inside a horizon volume. Practically all of 
them are eaten by the black holes. Mass density of 
s/Nm ~ 10^ monopoles can not be eliminated (under 
our conservative assumptions). However, the number of 
monopoles that the universe can tolerate is Nm ~ 10^. 
Thus, the mass density of the remnants does not violate 
observational constraints. 

Later, there can be further annihilation of monopole/ 
antimonopole pairs inside galaxies and clusters of galax- 
ies where these objects can clump; hence we shy away 
from extrapolating the freeze-out density in Eq.(l) (or 
Eq.(4)) all the way to the current epoch to estimate the 
mass dcinsity of remnants today and then applying the 
constraint of Eq(15). 

In conclusion, we have presented a possible solution to 
the cosmological monopole problem in which primordial 
black holes efficiently eliminate the magnetic monopoles 
from our universe. Under very conservative assumptions, 
black holes of mass 

Mbh < lO^gm (27) 

are dynamically capable of solving the cosmological 
monopole problem, while their required abundance does 
not violate any of the observational constraints. We did 
not assume any particular model for primordial black 
holes formation. Instead, we derived the bounds on 
abundance of black holes that can reasonably exist in 
our universe and are capable of removing monopoles. 
Black holes have a large cross-section for capturing non- 
rclativistic monopoles and the parameter space for solv- 
ing the monopole problem is not very restrictive. In a 
similar manner, black holes can in principle remove other 
unwanted massive relics from the early universe [18]. 
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